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Abstract 

We present the Fortran code SuSpect version 2.3, which calculates the Supersym¬ 
metric and Higgs particle spectrum in the Minimal Supersymmetric Standard Model 
(MSSM). The calculation can be performed in constrained models with universal 
boundary conditions at high scales such as the gravity (mSUGRA), anomaly (AMSB) 
or gauge (GMSB) mediated supersymmetry breaking models, but also in the non- 
universal MSSM case with R-parity and GP conservation. Gare has been taken to 
treat important features such as the renormalization group evolution of parameters 
between low and high energy scales, the consistent implementation of radiative elec- 
troweak symmetry breaking and the calculation of the physical masses of the Higgs 
bosons and supersymmetric particles taking into account the dominant radiative cor¬ 
rections. Some checks of important theoretical and experimental features, such as the 
absence of non desired minima, large fine-tuning in the electroweak symmetry break¬ 
ing condition, as well as agreement with precision measurements can be performed. 
The program is simple to use, self-contained and can easily be linked to other codes; 
it is rather fast and flexible, thus allowing scans of the parameter space with several 
possible options and choices for model assumptions and approximations. 


*The program with all relevant information can be downloaded from the web at the http site: 
www.lpta.univ-montp2.fr/~kneur/Suspect or obtained by sending an E-mail to one of the authors, 
abdelhak.dj ouadiOcern.ch, kneurSlpta.univ-montp2.fr, moultakaSlpta.univ-montp2.fr. 
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1. Introduction 


Supersymmetric theories (SUSY) [1], which provide an elegant way to stabilize the large 
hierarchy between the Grand Unihcation (GUT) and the electroweak scales and to cancel 
the qnadratic divergences of the radiative corrections to the Higgs boson masses, are by far 
the most stndied extensions of the Standard Model (SM). The most economical low-energy 
SUSY extension of the SM, the Minimal Snpersymmetric Standard Model (MSSM), which 
allows for a consistent nnihcation of the SM gange conplings and provides a natural solution 
of the Dark Matter problem, has been widely investigated; for reviews see Refs. [2-5]. As a 
corollary, the search for Snpersymmetric particles and for the extended Higgs spectrum has 
become the main goal of present and fnture high-energy colliders [6]. 

It is well-known that in the nnconstrained MSSM, it is a rather tedious task to deal 
with the basic parameters of the Lagrangian and to derive in an exhaustive manner their 
relationship with the physical parameters, i.e. the particle masses and conplings. This is 
mainly dne to the fact that in the MSSM, despite of the minimal gange gronp, minimal 
particle content, minimal couplings imposed by R-parity conservation and the minimal set 
of soft SUSY-breaking parameters, there are more than a hundred new parameters [7]. Even 
if one constrains the model to have a viable phenomenology [we will call later snch a model 
the phenomenological MSSM], assnming for instance no intergenerational mixing to avoid 
flavor changing nentral currents, no new source of GP violation, nniversality of hrst and 
second generation sfermions to cope with constraints from kaon physics, etc.., there are still 
more than 20 free parameters left. This large nnmber of inpnt enters in the evalnation of the 
masses of (P(30) SUSY particles and Higgs bosons as well as their complicated conplings, 
which involve several non-trivial aspects, such as the mixing between different states, the 
Majorana natnre of some particles, etc. The sitnation becomes particularly difficult if one 
aims at rather precise calculations and hence, attempts to include some rehnements such as 
higher order corrections, which for the calculation of a single parameter need the knowledge 
of a large part of, if not the whole, spectrum. 

Thus, the large number of free parameters in the nnconstrained or even phenomenological 
MSSM, makes a detailed phenomenological analysis of the spectra and the comparison with 
the outcome or expectation from experiment, a dannting task, if possible at all. Fortnnately, 
there are well motivated theoretical models where the soft SUSY-breaking parameters obey 
a nnmber of nniversal bonndary conditions at the high (GUT) scale, leading to only a 
handful of basic parameters. This is the case for instance of the minimal Snpergravity model 
(mSUGRA) [8], where it is assnmed that SUSY-breaking occurs in a hidden sector which 
commnnicates with the visible sector only through “flavor-blind” gravitational interactions. 
This leads to the simpler situation where the entire spectrnm of snperparticles and Higgs 
bosons is determined by the valnes of only Eve free parameters and makes comprehensive 
scans of the parameter space and detailed studies of the spectrum feasible. 
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However, there are also similarly constrained and highly predictive alternative SUSY- 
breaking models in the literature, such as anomaly mediated [9-11] or gauge mediated [12,13] 
SUSY-breaking models for instance, which should be investigated as well. We then have to 
trade a complicated situation where we have one model with many input parameters, with 
a not less complicated situation where we have many models with a small number of basic 
parameters. In addition, in these unihed models, the low-energy parameters are derived 
from the high-energy (GUT and/or possibly some intermediate scales) input parameters 
through Renormalization Group Equations (RGE) and they should also necessarily involve 
radiative electroweak symmetry breaking (EWSB), which sets additional constraints on the 
model. The implementation of the RG evolution and the EWSB mechanism poses numerous 
non-trivial technical problems if they have to be done in an accurate way, i.e. including 
higher order effects. This complication has to be added to the one from the calculation of 
the particle masses and couplings with radiative corrections (RG) which is still present. 

Therefore, to deal with the supersymmetric spectrum in all possible cases, one needs 
very sophisticated programs to encode all the information and, eventually, to pass it to 
other programs or Monte Garlo generators to simulate the physical properties of the new 
particles, decay branching ratios, production cross sections at various colliders, etc... These 
programs should have a high degree of flexibility in the choice of the model and/or the 
input parameters and an adequate level of approximation at different stages, for instance 
in the incorporation of the RGEs, the handling of the EWSB and the inclusion of radiative 
corrections to (super)particle masses, which in many cases can be very important. They 
should also be reliable, quite fast to allow for rapid comprehensive scans of the parameter 
space and simple enough to be linked with other programs. There are several public codes, 
in particular ISASUGRA [14], SOFTSUSY [15] and SPHENO [16], as well as a number of private 
codes, which deal with this problem. In this paper we present our program SuSpect. 

SuSpect, in its latest version 2.3 that we present here, is a Fortran code which calcu¬ 
lates the supersymmetric and Higgs particle spectrum in the constrained and unconstrained 
MSSMs. The acronym is an abbreviation of SusySpectrum and successive previous public 
versions of the code were available starting from 1997 and have been described in Ref. [17]. 
At the present stage, it deals with the “phenomenological MSSM” with 22 free parame¬ 
ters dehned either at a low or high energy scale, with the possibility of RG evolution to 
arbitrary scales, and the most studied constrained models, namely mSUGRA, AMSB and 
GMSB. Many “intermediate” models [e.g. constrained models but without unihcation of 
gaugino or scalar masses, etc..] are easily handled. The program includes the three major 
ingredients which should be incorporated in any algorithm for the constrained MSSMs: i) 
renormalization group evolution of parameters between a low energy scale [e.g. Mz and/or 
the EWSB scale] and a high-energy scale [e.g. Mqut], [17-19]; ii) consistent implementation 
of radiative electroweak symmetry breaking [loop corrections to the effective potential are 
included using the tadpole method] [20-23]; in) calculation of the physical (pole) masses 
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of the superparticles and Higgs bosons, including all relevant features such as the mixing 
between various states and the radiative corrections when important [24-34], 

The present version includes new options to read input files and write output hies in 
the recently adopted format of the SUSY Les Houches Accord (SLHA) for interfacing the 
spectrum generators with other computer codes [35]. The code contains three types of source 
hies: i) the main subroutine suspect2.f together with a separate routine twoloophiggs.f 
for the two-loop radiative corrections in the Higgs sector, ii) a separate calling routine hie 
suspect2_call.f, and Hi) two possible input hies suspect2.in (in the original SuSpect 
format) or suspect2_lha. in (in the SLHA format). Any choice and option is driven either 
from one of the two input hies [which is sufficient and convenient when dealing with a few 
model points] or alternatively directly from the suspect2_call. f hie, which also provides 
examples of calls for diherent model choices with all the necessary features [this option being 
useful when interfacing with other routines or when performing scans over the parameter 
space]. The program has several hags which allow to select the model to be studied and 
its input parameters, the level of accuracy of the algorithm [e.g. the iterations for the 
RGBs and the convergence of the EWSB], the level of approximation in the calculation 
of the various (s)particle masses [e.g. inclusion or not of RC], Besides the fact that it is 
hexible, the code is self-contained [the default version includes all routines needed for the 
calculation], rather fast [thus allowing large scans of the parameter space] and can be easily 
linked to other routines or Monte-Carlo generators [e.g. to calculate branching ratios, cross 
sections, relic densities]. All results, including comments when useful and some theoretical 
and experimental constraints, are found in the output hie suspect2.out (in the original 
SuSpect format by default) or alternatively in the output hie suspect2_lha.out, which are 
created at any run of the program. It is hoped that the code may be readily usable even 
without much prior knowledge of the MSSM. 

This “users’ manual” for the program, is organized as follows. In section 2, we briehy 
discuss the main ingredients of the unconstrained and phenomenological MSSMs as well as 
the constrained models mSUGRA, AMSB and GMSB, to set the notations and conventions 
used in the program. In section 3, we summarize the procedure for the calculation of the 
(s)particle spectrum: the soft SUSY-breaking terms [including the treatment of the input, 
the RG evolution and the implementation of EWSB], the physical particle masses [summa¬ 
rizing our conventions for the sfermion, gaugino and Higgs sectors]. We also discuss the 
theoretical [GGB, UFB, hne-tuning] and experimental [electroweak precision observables, 
muon g — 2, b ^ s'j branching fraction] constraints which can be imposed on the spectra, 
and how these are implemented in the code. In section 4, we summarize the basic practical 
facts about the program and discuss the content of the input and output hies with the possi¬ 
ble choices; we then make a brief comparison with other existing codes, discuss the interface 
with other programs and how the program is maintained on the web. A conclusion will be 
given in section 5. In the 
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2. The constrained and unconstrained MSSMs 

In this section, we summarize the basic assumptions which dehne the MSSM and the various 
constraints which can be imposed on it. This will also set the notations and conventions used 
in the program. We mainly focus on the unconstrained MSSM, the phenomenological MSSM 
with 22 free parameters, as well as on constrained models such as the minimal Supergrav¬ 
ity (mSUGRA), anomaly mediated (AMSB) and gauge mediated (GMSB) supersymmetry 
breaking models. 

2.1 The unconstrained MSSM 

The unconstrained MSSM is dehned usually by the following four basic assumptions [36]: 

(a) Minimal gauge group: the MSSM is based on the group SU(3)c x SU(2 )l x U(1)y, i-e. 
the SM symmetry. SUSY implies then that the spin-1 gauge bosons and their spin-1/2 
partners, the gauginos [bino B, winos lUi _3 and gluinos Gi-s], are in vector supermultiplets. 

(b) Minimal particle content: there are only three generations of spin-1/2 quarks and leptons 
as in the SM. The left- and right-handed chiral helds belong to chiral superhelds together 
with their spin-0 SUSY partners, the squarks and sleptons: Q,ur, dn, lJr. In addition, 
two chiral superhelds Hd, Hu with respective hypercharges —1 and -|-1 for the cancellation of 
chiral anomalies, are needed. Their scalar components, Hd and Hu, give separately masses 
to the isospin -1-1/2 and —1/2 fermions and lead to hve Higgs particles: two GP-even h,H 
bosons, a pseudoscalar A boson and two charged H^ bosons. Their spin-1/2 superpartners, 
the higgsinos, will mix with the winos and the bino, to give the “ino” mass eigenstates: the 
two charginos xf 2 the four neutralinos Xi 2 3 4 - 

(c) Minimal Yukawa interactions and R-parity conservation: to enforce lepton and baryon 
number conservation, a discrete and multiplicative symmetry called R-parity is imposed. It 
is dehned by Rp = (_i^2s+3B+l^ where L and B are the lepton and baryon numbers and s 
is the spin quantum number. The R-parity quantum numbers are then Rp = +1 for the 
ordinary particles [fermions, gauge and Higgs bosons], and Rp = —1 for their supersymmetric 
partners. In practice, the conservation of R-parity has important consequences: the SUSY 
particles are always produced in pairs, in their decay products there is always an odd number 
of SUSY particles, and the lightest SUSY particle (LSP) is absolutely stable. 

The three conditions listed above are sufficient to completely determine a globally super- 
symmetric Lagrangian. The kinetic part of the Lagrangian is obtained by generalizing the 
notion of covariant derivative to the SUSY case. The most general superpotential, compat¬ 
ible with gauge invariance, renormalizability and R-parity conservation is written as: 

^= 5 ^ -YtjUmHu.Qj + Yf^dR,Hd.Qj + YyiR,Hu.Lj + fiHu.Hd ( 1 ) 

i,j=gen 

The product between SU(2 )l doublets reads H.Q = CabH'^Q^ where a,b are SU(2 )l indices 
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and ei 2 = 1 = —£ 21 , and denote the Yukawa couplings among generations. The first 

three terms in the previous expression are nothing else but a superspace generalization of 
the Yukawa interaction in the SM, while the last term is a globally supersymmetric Higgs 
mass term. The supersymmetric part of the tree-level potential Vtree is the sum of the 
so-called F- and D-terms [37], where the F-terms come from the superpotential through 
derivatives with respect to all scalar helds 0a, Vp = with fF“ = dW/d<pa, and the 

D-terms corresponding to the U(1)y, SU(2)l, and SU(3)c gauge symmetries are given by 
Vd = \ with T* and gi being the generators and the coupling constants 

of the corresponding gauge groups. 

(d) Minimal set of soft SUSY-breaking terms: to break Supersymmetry, while preventing the 
reappearance of the quadratic divergences [soft breaking], one adds to the supersymmetric 
Lagrangian a set of terms which explicitly but softly break SUSY [38]: 

• Mass terms for the gluinos, winos and binos: 

3 8 

MiBB + M 2 ^W^Wa + M^^G^Ga + h.c. 

a=l a=l 

• Mass terms for the scalar fermions: 

-^sfermions ^ ^ ^^^Q^QiQi T T ^l^Ri\ (3) 

i=gen 

• Mass and bilinear terms for the Higgs bosons: 

“ -^Higgs = + mlHlH, + + h.c.) (4) 

• Trilinear couplings between sfermions and Higgs bosons 

-A,I 1 .= Y + + h.c.] (5) 

i,j=gen 

The soft SUSY-breaking scalar potential is the sum of the three last terms: 

hsoft -^sfermions -^Higgs -^tril. (6) 

Up to now, no constraint is applied to this Lagrangian, although for generic values of the pa¬ 
rameters, it might lead to severe phenomenological problems such as flavor changing neutral 
currents [FCNC] and unacceptable amount of additional CP-violation [39], color and/or 
charge breaking minima [40], an incorrect value of the Z boson mass, etc... The MSSM 
dehned by the four hypotheses {a)-{d) above, will be called the unconstrained MSSM. 
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2.2 The “phenomenological” MSSM 

In the unconstrained MSSM, and in the general case where one allows for intergenerational 
mixing and complex phases, the soft SUSY breaking terms will introduce a huge number 
(105) of unknown parameters, in addition to the 19 parameters of the SM [7]. This large 
number of free parameters makes any phenomenological analysis in the general MSSM very 
complicated as mentioned previously. In addition, many “generic” sets of these parameters 
are excluded by the severe phenomenological constraints discussed above. A phenomenolog¬ 
ically viable MSSM can be dehned by making the following three assumptions: [i) All the 
soft SUSY-breaking parameters are real and therefore there is no new source of CP-violation 
generated, in addition to the one from the CKM matrix, [ii) The matrices for the sfermion 
masses and for the trilinear couplings are all diagonal, implying the absence of FCNCs at the 
tree-level. {Hi) First and second sfermion generation universality at low energy to cope with 
the severe constraints from mixing, etc [this is also motivated by the fact that one 

can neglect for simplicity all the masses of the hrst and second generation fermions which 
are too small to have any effect on the running of the SUSY-breaking parameters]. 

Making these three assumptions will lead to 22 input parameters only: 
tan/5: the ratio of the vevs of the two-Higgs doublet helds. 

the Higgs mass parameters squared. 

Ml, M 2 , M 3 : the bino, wino and gluino mass parameters. 

the first/second generation sfermion mass parameters. 
ttIq, 171^^,1711^, mi, nirj^: the third generation sfermion mass parameters. 
Au,Ad,Ae. the first/second generation trilinear couplings. 

At, Ab, Ar'. the third generation trilinear couplings. 

Two remarks can be made at this stage: {i) The Higgs-higgsino (supersymmetric) mass 
parameter |/u| (up to a sign) and the soft SUSY-breaking bilinear Higgs term B are deter¬ 
mined, given the above parameters, through the electroweak symmetry breaking conditions 
as will be discussed later. Alternatively, one can trade the values of m‘\j^ and with the 
“more physical” pseudoscalar Higgs boson mass Ma and parameter fi [such an alternative 
choice is explicitly possible in SuSpect by appropriate setting of the input parameters], {ii) 
Since the trilinear sfermion couplings will be always multiplied by the fermion masses, they 
are important only in the case of the third generation. However, there are a few (low scale) 
situations, such as the muon {g — 2) and the neutralino-nucleon scattering for direct Dark 
Matter searches, where they will play a role. We thus consider them as input. 

Such a model, with this relatively moderate number of parameters [especially that, in 
general, only a small subset contributes dominantly (i.e. at tree-level) when one looks at a 
given sector of the model] has much more predictability and is much easier to investigate 
phenomenologically, compared to the unconstrained MSSM. We will refer to this 22 free 
input parameters model as the “phenomenological” MSSM or pMSSM [4]. 
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2.3 The mSUGRA model 


Almost all problems of the general or unconstrained MSSM are solved at once if the soft 
SUSY-breaking parameters obey a set of universal boundary conditions at the GUT scale. 
If one takes these parameters to be real, this solves all potential problems with CP violation 
as well. The underlying assumption is that SUSY-breaking occurs in a hidden sector which 
communicates with the visible sector only through gravitational-strength interactions, as 
specihed by Supergravity. Universal soft breaking terms then emerge if these Supergravity 
interactions are “flavor-blind” [like ordinary gravitational interactions]. This is assumed to 
be the case in the constrained MSSM or minimal Supergravity (mSUGRA) model [ 8 ]. 

Besides the unihcation of the gauge coupling constants ( 71 , 2,3 of the U(l), SU(2) and SU(3) 
groups, which is well verihed given the experimental results from LEPl [41] and which can 
be viewed as hxing the Grand Unihcation scale Mqut 2 ■ 10^® GeV [42], the unihcation 
conditions in mSUGRA, are as follows: 

- Unihcation of the gaugino [bino, wino and gluino] masses: 

Afi(MGUT) = Af 2 (MGUT) = Af3(MGUT) = ''^ 1/2 (7) 

- Universal scalar [i.e. sfermion and Higgs boson] masses [i is the generation index]: 

^Qi(^GUT) = Mu^,(Mgut) = Mj^^(Mgut) = M^,(Mgut) = M|-^,(Mgut) 

= MhSMq\]-£) = = mo (8) 

~ Universal trilinear couplings: 

^p(^gut) = = AqSij (9) 

Besides the three parameters mi/ 2 , uiq and Aq, the supersymmetric sector is described at 
the GUT scale by the bilinear coupling B and the supersymmetric Higgs(ino) mass parameter 
jj,. However, one has to require that EWSB takes place at some low energy scale. This results 
in two necessary minimization conditions of the two-Higgs doublet scalar potential which, 
at the tree-level, has the form [to have a more precise description, one-loop corrections to 
the scalar potential have to be included, as will be discussed later]: 

Uffiggs = m\HlHd + mlHlHu + + h.c.) 



where we have used the usual short-hand notation: ml = + /i^, 

m 3 = Bfi and the SU(2) invariant product of the two doublets ■ 02 = 0 i 02 ~ 0 i 02 - The 
two minimization equations dVmggs/dH^ = dVmggs/dH^ = 0 can be solved for and Bjj,-. 

^ tan 2/5(m^^ tan f3 — cot 0) — M| 

R/i = ^ sin 20 m]j^ + + 2fi^ 
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( 11 ) 



(which are of course valid as well in the more general unconstrained MSSM). Here, M| = 
{di + S'!) ■ ('^« + /4 and tan (3 = Vu/vd is dehned in terms of the vacuum expectation values 

of the two neutral Higgs helds. Consistent EWSB is only possible if eq. ca gives a positive 
value of /i^. The sign of /x is not determined. Therefore, in this model, one is left with only 
four continuous free parameters, and an unknown sign: 

tan/3 , mi /2 , mo , Hq , sign(/i). (12) 

All the soft SUSY breaking parameters at the weak scale are then obtained through Renor¬ 
malization Group Equations. 


2.4 The AMSB model 


In mSUGRA, Supersymmetry is broken in a hidden sector and the breaking is transmitted 
to the visible sector by gravitational interactions. In Anomaly Mediated Supersymmetry 
Breaking models, the SUSY-breaking occurs also in a hidden sector, but it is transmitted 
to the visible sector by the super-Weyl anomaly [9,10]. The gaugino, scalar masses and 
trilinear couplings are then simply related to the scale dependence of the gauge and matter 
kinetic functions. This leads to soft SUSY-breaking scalar masses for the hrst two generation 
sfermions that are almost diagonal [when the small Yukawa couplings are neglected]. 

In terms of the gravitino mass m 3/2 [which is much larger than the gaugino and squark 
masses, a cosmologically appealing feature], the f3 functions for the gauge and Yukawa cou¬ 
plings Qa and Yj, and the anomalous dimensions y* of the chiral superhelds, the soft SUSY 
breaking terms are given by: 


Ma 

A. 


9a 

Yi 


W3/2 


■m3/2 


-- ( E, 


d-fi 

dga 


A 






^ 1/2 


(13) 


These equations are RG invariant and thus valid at any scale and make the model highly 
predictive. The additional parameters, /i^ and B are obtained as usual by requiring the 
correct breaking of the electroweak symmetry. One then has, in principle, only three input 
parameters 771 ^/ 2 ,tanj3 and sign(p). However, this rather simple picture is spoiled by the 
fact that the anomaly mediated contribution to the slepton scalar masses squared is negative 
and the sleptons are in general tachyonic. This problem can be cured by adding a positive 
non-anomaly mediated contribution to the soft masses. The simplest phenomenological way 
of parameterizing the non-anomaly contribution is to add a common mass parameter mo 
at the GUT scale, which would be then an additional input parameter to all the (squared) 
scalar masses. However in the general case, the non-anomaly mediated contribution might 
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be different for different scalar masses and depend on the specific model which has been 
chosen. One should then write a general non-anomalous contribution at the GUT scale for 
each scalar mass squared: 

= cs,ml - ^ + D terms. (14) 

where the coefficients c^. depend on the considered model. 

A few examples of models with different non-anomalous contributions are: 

- The minimal anomaly mediated supersymmetry breaking model with a universal mo [43]: 

Cq Cl Cgjj CHy Ch^ 1 

- The gaugino assisted AMSB model where one assumes that gauge and gaugino helds reside 
in the bulk of an extra dimension [44]: 

cq = 21/10, = 8/5, = 7/5, cl = 9/10, Ce = 3/5, ch^ = 9/10ch^ (16) 

- Models where an extra U(l) factor is added; a particular scenario is interesting phenomeno¬ 
logically since it leads to a light top squark [45]: 

Cq 3, Cyj^ CcIr 1) Cl Cg 1, Cjj^ Ch^ 2 (l^^) 

A simple way to account for all the different models is to add to the three continuous 
and one discrete original basic parameters, the set of coefficients c^. as input to specify, and 
therefore one would have the set of input parameters: 

mo , mo /2 , tan/3 , sign(p) and C 5 . (18) 

This is the approach that we will follow in the program. 

2.5 The GMSB model 

In Gauge Mediated Supersymmetry Breaking models, SUSY-breaking is transmitted to the 
MSSM helds via the SM gauge interactions. In the original scenario [46], the model consists 
of three distinct sectors: a secluded sector where SUSY is broken, a “messenger” sector con¬ 
taining a singlet held and messenger helds with SU(3)c x SU(2)l x U(1)y quantum numbers, 
and a sector containing the helds of the MSSM. Another possibility, the so-called “direct 
gauge mediation” [47] has only two sectors: one which is responsible for the SUSY breaking 
and contains the messenger helds, and another sector consisting of the MSSM helds. In 
both cases, the soft SUSY-breaking masses for the gauginos and squared masses for the 
sfermions arise, respectively, from one-loop and two-loop diagrams involving the exchange 
of the messenger helds, while the trilinear Higgs-sfermion-sfermion couplings can be taken 
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to be negligibly small at the messenger scale since they are [and not their square as for the 
sfermion masses] generated by two-loop gauge interactions. This allows an automatic and 
natural suppression of FCNC and CP-violation; for a review see, Ref. [12]. 


In the GMSB models that we will consider, the source of SUSY breaking is parameterized 
by an SU(3)c x SU(2)l x U(1)y gauge-singlet chiral superheld S whose scalar and auxiliary 
components acquire vacuum expectation values denoted by S and Fs, respectively. We 
assume rig pairs of q, q quark-like [resp. Uj- pairs of /, I lepton-like] messenger superhelds 
transforming as (3,1, —I), (3,1, |) [resp. (1, 2, |), (1, 2, —|)] under SU(3)c x SU(2 )l x U(1)y 
and coupled to S through a superpotential of the form XSqq + XSU. Soft SUSY-breaking 
parameters are then generated at the messenger scale = ^S, 






.^/(Vmes) 



(19) 

( 20 ) 
( 21 ) 


where A = Fs/S, G = U(l), SU(2), SU(3), m labels the messengers and s runs over the 
Higgs doublets as well as the left-handed doublets and right-handed singlets of squarks and 
sleptons. The one- and two-loop functions g and / are given by [Li 2 is the Spence function]: 


9{x) 

fix) 



1 + X 


+ x) log(l + a:) + (1 — x) log(l — a:)] 

_log(l + .)-2Li,(^)+lLi,(^) 


+ {x —x) 


( 22 ) 


Defining the Dynkin index by 


Tr (T^rJ) 


2 


^ab 


(23) 


for non-abelian groups, and = (6/5)Y^ where Y = Qem — one has (see ea. (IT!H) l 


= ^(2^4+ 3n;) 

= ni 

(24) 


With the Casimir invariant given by 

= CS"'"’! = ^^^1 (25) 

for the N of SU{N), and = (3/5) Y^, one hnds for 


AfC{s) 


E 

m,G 


Of(^(A/nies) 

Air 


V|(m)Gf(.) 
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(see eq. (EDD) the following values: 


UC{Q) 

uc{ij) 

UC{D) 

UC{L) 

UC{E) 

J^C{Hu) 


IGtt^ 

1 

IGtt^ 

1 

IGtt^ 

1 

167r^ 

1 

167r^ 


- Tij \ 9 ^Tlj 9 

—h H-H-+ 

100 150^ ^ 4 ^ 


,4nr 8 ng, 9 4ng 9 

, Tlj 2 ??-^ \ 2 ‘^^(5 2 

^9^ I 2 

^100+ 50^^+ 4 


4n, 


Q 2 

ao 


_i_ ^^9 ^^.2 

^ 25 25 ’ ^ 


UC{Hd) = AfC{L) 


(26) 


The freedom in choosing independently the number of nq and rif messengers allows to study 
various model conhgurations: for instance when the messengers are assumed to form com¬ 
plete representations of some grand unihcation group (e.g. 5 -|- 5 of SU(5)) where riq = rif, or 
when they transform under larger unihcation group factors with some extra discrete symme¬ 
tries where typically riq 7 ^ nq [48]. [When nq = Uf = 1 one retrieves the minimal model [46]. 
In this case the gaugino masses have the same relative values as if they were unihed at 
Mgut despite the fact the boundary conditions are set at Mmes and that scalar masses are 
havor independent. Furthermore when A/Mmes "C 1, one has f{x) ~ g{x) ~ 1 .] In addition, 
some constraints are in general needed in order to have a viable spectrum, for instance: 
A/Mmes < 1 to avoid negative mass squared for bosonic members of the messenger scale and 
A/Mmes ^ 0.9 to avoid too much hne-tuning in EWSB. Note also that nq > nf improves the 
hne-tuning issue [49]. 

Once the boundary conditions are set at the low energy parameters are obtained via 

the usual RGBs and the proper breaking of the EW symmetry is required. 

Therefore, in the GMSB model that we are considering, there are six input parameters 


tan/? , sign(/i) , M^es , A , nq , nq (27) 

In addition, one has to include as input the mass of the gravitino G which, in this case 
is the lightest SUSY particle. This mass, = F/{\/3Mp) with Mp the reduced Planck 
mass, will depend on an additional free parameter F which parameterizes the scale of the 
full SUSY breaking and whose typical size is of 0{Fs) in direct mediation and much larger 
in secluded mediation. The choice of this parameter, which plays a role only for the lifetime 
of the next-to-lightest SUSY particle, is left to the user. 


2.6 Non—universality models 

mSUGRA, AMSB and GMSB are well dehned models of which the possible phenomenolog¬ 
ical consequences and experimental signatures have been widely studied in the literature. 
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However, in the absence of a truly fundamental description of SUSY-breaking, none of these 
models should be considered as the compelling model. Unknown physics at the Planck 
and/or GUT scales renders some of the basic assumptions inherent to the above scenarii 
rather uncertain: deviations from mSUGRA universality conditions at the high scale are 
generally expected in non-minimal Supergravity [50], or in superstring [51] settings while de¬ 
viations from the initial AMSB scenario can be expected, due to extra anomaly contributions 
coming from the underlying fundamental theory [10] and deviations from minimal GMSB 
can occur due for instance to the breaking of symmetries which protected direct coupling 
between messenger and MSSM matter helds. 

To be on the safe side from the experimental point of view, it is therefore wiser to allow for 
a departure from these models, and to study the phenomenological implications of relaxing 
some dehning assumptions. However, it is often desirable to limit the number of extra 
free parameters, in order to retain a reasonable amount of predictability when attempting 
detailed investigations of possible signals of SUSY. Therefore, it is more interesting to relax 
only one [or a few] assumption[s] at a time and study the phenomenological implications. 
Of course, since there are many possible directions, this would lead to several intermediate 
MSSMs between these constrained models and the phenomenological MSSM with 22 free 
parameters discussed in section 2.2. 

Taking the most studied mSUGRA as a reference model, we briefly discuss here some 
model cases where the universality conditions are naturally violated: 

i) non unihcation of the soft SUSY-breaking gaugino mass terms: 

M,{Mu) + M^iMu) + (28) 

This occurs for instance in Superstring motivated models in which the SUSY breaking is 
moduli dominated such as in the O-I and 0-H models [52], or in extra dimensional SUSY- 
GUT models in which the additional dimensions lead to the breaking of the large gauge 
symmetry and/or supersymmetry or in SUSY models where the breaking occurs through a 
non SU(5) singlet F term; see Ref. [53] for phenomenology oriented discussions. 

ii) mSUGRA with non-unihcation of the two first and third generation scalar masses: 

^OQ = ^0 L---7^ = ^01 ■ ■ ■ (29) 

This occurs in models where the soft SUSY-breaking scalar masses at the GUT scale are 
influenced by the fermion Yukawa couplings. This is the case for the so-called inverted mass 
hierarchy models [54] where the scalar mass terms of the hrst two generations can be very 
heavy (9(10 TeV), while those of the third generation sfermions and the Higgs bosons are 
rather light, solving thus the SUSY flavor and GP problems, which are related to the first 
two generations, while still satisfying naturalness constraints. 

Hi) mSUGRA-like models, but with non-universality of the sfermion and Higgs boson 
scalar masses [keeping a priori universal gaugino masses and trilinear soft breaking cou¬ 
plings] . In a model-independent context, one can in principle consider any non-universality 
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configuration which does not violate important constraints from FCNC, etc... Interestingly, 
such configurations occur naturally in models where the grand unification group has a rank 
strictly greater than 4 [55]. Starting from universal soft scalar masses, non-universality ef¬ 
fects occur at intermediate scales, related to the spontaneous symmetry breaking down to 
the standard model gauge group, via D-term contributions corresponding to the extra U(l) 
symmetries [56]. For instance, in SO(IO) SUSY GUT models with universal boundary con¬ 
ditions, and assuming that the GUT group breaks directly to the SM group at the GUT 
scale, one is lead to patterns of the form [56,57], 


= ^eu = ^un + = mi^ Mh^ ^ Mh^ (30) 

where there are actually three free mass parameters, two of which can be taken as the 
initial common soft masses for the sfermion and Higgs SO(IO) multiplets, the third being 
the D-term contribution associated to the reduction of the group rank. 

In the present SuSpect version such specific SUSY GUT models are not yet available^, 
but one can still have a reasonable estimate of the non-universality effects by using the 
pMSSM option. Alternatively, one can also use the SuSpect option which allows to discon¬ 
nect the Higgs sector from the sfermionic by introducing two additional input parameters: 
the pseudoscalar Higgs boson mass Ma and the higgsino mass parameter fi [which have a 
more direct “physical” interpretation than the scalar mass terms Mh„, This allows to 

perform more general phenomenological or experimental analyses; c.f. some LEP and LHG 
Higgs analyses [58] or some recent Dark Matter studies [59]. 

iv) Partially unified models where one relaxes one or a few parameters to fit some collider 
zoo event or to analyze a phenomenological situation which introduces new features. This 
is the case, for instance, for the light top squark scenario which can be set by hand to 
discuss some theoretical issues [such as baryogenesis in the MSSM [60] for instance] or 
phenomenological situations [such as new decay or production modes of top squarks [61] 
for instance] . 

An easy and practical way to implement these various non-unified or partially unified 
scenarii, is to allow for the possibility of choosing all the soft SUSY-breaking parameters 
listed above for the phenomenological MSSM of section 2.2 [the 22 parameters except for 
tan/?] at the high-energy or GUT scale, with the boundary conditions set by hand and 
chosen at will. One can even chose the scale at which the boundary conditions are set to 
account for intermediate scales. If this scale is the electroweak symmetry breaking scale, 
then we have simply the MSSM with the soft SUSY-breaking parameters defined at the low 
energy scale, i.e. the phenomenological MSSM. All these options are provided by our code. 

^These would require modified RGBs at least for some sectors, inclusion of the effects of new heavy fields 
(right-handed (s)neutrinos for instance), of threshold effects corresponding to the various scales of successive 
symmetry breaking down to the SM, etc... 
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3. The Particle Spectrum Calculation with Suspect 


In this section, we discuss our procedure for calculating the SUSY and Higgs particle spec¬ 
trum. We will take as example the sophisticated cases of the constrained MSSMs with 
universal boundary conditions at the high scale, mSUGRA AMSB and GMSB, where all 
ingredients included in the SuSpect algorithm are present: RGBs, radiative EWSB and cal¬ 
culation of the physical particle masses. We hrst describe the general algorithm, then discuss 
the calculation of the soft SUSY-breaking terms, the determination of the particle masses, 
the various theoretical and phenomenological tests that we impose on the model parameters. 

3.1 General algorithm 

As mentioned previously, there are three main steps for the calculation of the supersymmetric 
particle spectrum in constrained MSSMs, in addition to the choice of the input parameters 
and the check of the particle spectrum: 

i) Renormalization group evolution (RGE) of parameters [17-19], back and forth between 
the low energy scales, such as Mz and the electroweak symmetry breaking scale, and the 
high-energy scale, such as the GUT scale or the messenger scale in GMSB models. This RGE 
is performed for the SM gauge and Yukawa couplings and for the soft SUSY-breaking terms 
(scalar and gaugino masses, bilinear and trilinear couplings and tan j3) and /i. This procedure 
has to be iterated in order to include SUSY threshold effects or radiative corrections due to 
Higgs and SUSY particles. 

ii) The consistent implementation of (radiative) electroweak symmetry breaking [20-23] 
and the calculation of B and |/i| from the effective scalar potential at one-loop level (plus 
the leading two-loop contributions). Here, we use the tadpole method to include the loop 
corrections [22]. The procedure has to be iterated until a convergent value for these two 
parameters is obtained. [In the hrst step, the values of and the electroweak symmetry 
breaking scale are guessed by using the tree-level potential since no sparticle or Higgs mass 
has been calculated yet.] 

in) Galculation of the pole masses of the Higgs bosons and the SUSY particles, including 
the mixing between the current states and the radiative corrections when they are important 
[24-31]. In this context, we will follow mainly the content and notations of [25], to which 
we will refer as PBMZ. The latter provides most of the necessary contributions at the one- 
loop level for the Higgs and sparticle masses, while the leading two-loop corrections to the 
(neutral) Higgs masses, calculated in Refs. [30-33], are implemented in the form of a separate 
subroutine twoloophiggs. f (provived by P. Slavich). 

The general algorithm is depicted in Figure 1, and we will discuss the various steps in some 
detail in the following subsections. 
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Figure 1: Iterative algorithm for the ealeulation of the SUSY partiele speetrum in SuSpect 
from the ehoice of input (first step) to the eheck of the speetrum (last step). The steps 
are detailed in the various subseetions. The EWSB iteration [ealeulationally fast] on p, 
is performed until \pi — Pi-i\ < e\pi\ (with e ~ while the RG/RC “long” iteration 

[ealeulationally longer] is performed until some (user speeified) aeeuraey is reaehed 
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3.2 Calculation of the MSSM parameters at the low scale 

Prior to the calculation of all the relevant terms entering the physical mass spectrum calcula¬ 
tion at a low (EWSB) scale, such as typically the soft SUSY-breaking terms in constrained 
models with boundary conditions at the unification scale, the first important stage is to 
define a choice of low-energy input parameters and how they are extracted from the present 
experimental data. This is the purpose of the next subsection. 


3.2.1 Choice and treatment of the SM inputs 


The low-energy (weak) scale boundary conditions set the gauge and Yukawa couplings by 
matching the running MSSM parameters to the experimental data at a chosen renormaliza¬ 
tion scale, usually Q = Mz- This step is quite involved, requiring to subtract the radiative 
corrections from the experimental data in order to arrive at the DR-renormalized MSSM 
parameters. Using the formulae of mainly Ref. [25] and references therein, the MSSM DR 
gauge couplings gi ,g 2 , gs and the electroweak parameter v can be computed at Q = Mz 
from a set of four experimental input parameters. 

The latter can be chosen as: Gp, the Fermi constant determined from the muon decay; 
Mz, the pole mass of the Z boson; aem{Mz)^^, the five-flavour SM electromagnetic coupling 
at the scale Mz in the MS scheme; as{Mz)^^, the five-flavour SM strong coupling at the 
scale Mz in the MS scheme. Then, the running couplings g 2 , gi are connected to the running 
Z-boson mass mz by the relation: 

ml = Ml + KelllPMl) = j (s? + gl) vl (31) 

where is the transverse part of the Z boson self-energy computed at a squared 

external momentum equal to the squared pole Z boson mass. [Note that according to this 
input choice, the W boson mass and the parameter sin^ 6 w are not free parameters but can 
be derived when needed, e.g. in the DR scheme, from the above dehned input parameters 
using appropriate relations including radiative corrections]. 

The next MSSM input parameter is tan/? = Vu/vd, the ratio of the vacuum expectation 
values (VEVs) of the two MSSM neutral Higgs helds, and H^, dehned also at the scale 
Q = Mz, with vl + vl = v\ 

Then the Yukawa couplings {u = u, c, t) for the up-type quarks. Yd {d = d, s, h) for 
the down-type quarks and = e, fi, r) for the leptons are determined in the DR scheme 
from the corresponding running fermion masses as 


^/2m^ ^ _ V2md ^ _ y/2me 

n sin /? ’ ^ V cos /? ’ ^ v cos /? ’ 


(32) 


The DR running fermion masses (see e.g. Refs. [62-65]) m/ (with / = u,d,i) in eq. 
can be derived at the one-loop level from the corresponding pole masses Mj through the 
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relation 


mf = Mf + Ef{Mf) (33) 

where S/(M/) is the one-loop fermion self-energy computed at an external momentum equal 
to the pole mass. In the case of the top quark, the self-energy includes also the leading two- 
loop standard QCD corrections [63] 

where the precise expression for depends on the renormalization scheme in 

which the parameters entering the one-loop self-energy are expressed. Concerning the 
bottom quark, we follow now the SUSY Les Houches Accord [35] which prescribes to take 
as input the SM running mass in the MS scheme, In addition, a “resummation” 

procedure is required (see e.g. Ref. [66]) in order to properly take into account the large QCD 
corrections, as well as the tan/3-enhanced SUSY corrections [67], to the relation between 
the input bottom mass and the corresponding MSSM, DR Yukawa coupling. We extract the 
latter, via eq. (ISS), from the MSSM, DR bottom mass fhb, defined at the scale Q = Mz by 
the following matching condition: 

mb = '^b{Mz)^sM — I _ 

where nib = ’^b(^z)sM ^R bottom mass, obtained by evolving mb{mb)^^ up 

to the scale Q = Mz with the appropriate RGE, in order to resum the QCD corrections, 
and then converting it to the DR scheme; Ab = Tjb{mb)/nib accounts for the remaining non- 
gluonic corrections, some of which are enhanced by a factor tan/3. It has been shown [66] 
that defining the running MSSM bottom mass as in eq. (jnsi) guarantees that the large 
threshold corrections of 0{asianij3)'^ are included in fhb to all orders in the perturbative 
expansion. Concerning the r lepton, the only tan /3-enhanced corrections to be included are 
those controlled by the electroweak gauge couplings, from chargino-slepton loops. 

Concretely the (present) reference values of the SM input parameters at the weak scale 
are hxed as follow. We take for the electroweak and strong parameters [41]: 

Gf = 1.16639 10-® GeV^ Mz = 91.1876 GeV, 

af^{Mzf^ = 127.934, as{Mzf^ = 0.1172, (36) 

and for the third-generation fermion masses the values [41,68]: 

Mt = 178.0 GeV, mb{mbf^ = 4.25 GeV, M, = 1.777 GeV. (37) 

Of course, for maximal flexibility, all those SM input default values in eqs. ()36l37j) can be 
changed at will in the input file. 
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Once the Supersymmetric particle spectrum has been obtained [see below], we include 
all the important SUSY radiative corrections to the gauge and Yukawa couplings, via their 
above dehning relations from input parameters, e.g. eqs. mm or (jnSD, where the SUSY 
particle masses enter the expression of e.g. and in eq. (IH^ . This neces¬ 

sarily implies an iterative procedure, since the values of the SUSY particle spectrum depend, 
among other thing, on the precise values of the gauge and Yukawa couplings via the RGE 
typically. The iteration is to be performed until a sufficiently stable hnal SUSY spectrum 
is obtained [and where the very Erst iteration is done with only SM radiative corrections 
included in eqs. mm . before SUSY particle masses are dehned, see the general algorithm 
overview in Fig. 1]. In the case of the Yukawa couplings, we include all relevant SUSY 
corrections to the third generation fermion masses. For the bottom quark (r lepton) mass, 
we include the SUSY-QCD and stop-chargino (sneutrino-chargino) one-loop corrections at 
zero-momentum transfer [66] which, according to PBMZ, is an extremely good approxima¬ 
tion. These corrections to the b and r masses are enhanced by terms oc ptan/3 and can 
be rather large. For the top quark, the inclusion of only the leading corrections at zero 
momentum transfer is not an accurate approximation, and we include the full one-loop 
SUSY-QCD [i.e stop and gluino loops] and electroweak [i.e. with gauge, Higgs boson and 
chargino/neutralino exchange] corrections d la PBMZ [25]. 

3.2.2 Renormalization Group Evolution 

All gauge and (third generation) Yukawa couplings are then evolved up to the GUT scale 
using the two-loop MSSM RGBs [19, 20] in the DR scheme, with the contribution of all 
the MSSM particles in the relevant beta-functions. The GUT scale, Mgut — 2 ■ 10^® GeV 
can be either Exed by hand or, by appropriate user’s choice in the input Ele, calculated 
consistently to be the scale at which the electroweak gauge coupling constants [with the 
adequate normalization] unify, gi = g 2 ■ ^"3/5. In contrast, we do not enforce exact g 2 = gs 
uniEcation at the GUT scale and assume that the small discrepancy, of at most a few percent, 
is accounted for by unknown GUT-scale threshold corrections [69]. 

One can then chose the parameter tan/?, given at the scale Mz, the sign of the /i pa¬ 
rameter and, depending on the chosen model, the high energy and the low energy input. 
For instance, one can set the high-energy scale Enigh, which in mSUGRA or AMSB can be 
either forced to be Mqut [the scale at which gi and g 2 unify] or chosen at will [any particular 
intermediate scale between Mz and Mgut can be allowed in general and in the case of the 
GMSB model this scale corresponds to the messenger scale Mmes]- Similarly the low energy 
scale Flow, where the RGFs start or end may be chosen [it is in general taken to be Mz or 
the EWSB scale to be discussed later]. The additional input in the various models are: 

• mSUGRA : the universal trilinear coupling Aq, the common scalar mass mo and the 
common gaugino mass mi/ 2 , all deEned at the scale Mqut- 
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• AMSB : the common scalar mass mo, the gravitino mass m 3/2 and the set of coefficients 
csi for the non-anomalous contributions, to be as general as possible. 

• GMSB : the scale A, the messenger scale Mmes which corresponds to Emgh, as well as 
the numbers of messengers Ug and n;. 

• pMSSM with boundary conditions: the various soft SUSY-breaking parameters listed 

in section 2.2 [21 parameters in total, in addition to tan/?] defined at the scale A’mgh- 
These input can also be chosen at will at the low-energy scale i?Low which has also 
to be provided as a necessary input. For this pMSSM input choice we stress that all 
the (scale-dependent) soft SUSY-breaking parameters are implicitly understood to be 
dehned at the given arbitrary low scale Ulow In particular, the code will perform 
RGE consistently between the scale Q = Mz and Q = Ulow, to match e.g. with the 
SM input dehned at Q = mz, as well as the RGE between Q = Elow and the EWSB 
scale Mewsb, while the RGE from a high scale Enigh down to Mewsb are not relevant 
in this case and are thus switched off. Note also that here, a very convenient option 
is provided which allows to trade the input parameters and with the more 
“physical” parameters M 4 (pole mass) and [again in such a way that EWSB 

is consistently realized, with a warning hag whenever it is not the case]. 

Given these boundary conditions, all the soft SUSY breaking parameters and couplings 
are evolved down to the EWSB and Mz scales, using either the one- or (preferably) full 
two-loop RGE options^. Our default choice for the EWSB scale is the geometric mean of 
the two top squark running masses (in the DR scheme), 

^^EWSB = (38) 

which minimizes the scale dependence of the one-loop ehective potential [23] discussed below 
[at hrst iteration where the stop masses have not yet been calculated, we use the geometric 
mean of the soft SUSY-breaking stop masses instead as a hrst guess]. Note, however, that 
any other arbitrary values of the EWSB scale can be chosen easily by an appropriate input 
setting. Since tan/? is dehned at Mz, the vevs have to be evolved down from Mewsb to Mz- 

Once the SUSY spectrum is calculated [see below], the heavy (s)particles are taken to 
contribute to reevaluate the gauge and Yukawa couplings at the scale Q = Mz, as discussed 
in the previous subsection, and the necessary iterative procedure thus also includes the 
RGEs. 

^The full two-loop RGEs in the MSSM [70], have been implemented for the slepton and quark masses in 
the latest version of SuSpect . In the previous versions, the two-loop RGEs were included for all terms except 
for the latter soft SUSY-breaking masses. These higher order effects can have a non-negligible impact [95] 
in some regions of the MSSM parameter space. 
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3.2.3 Electroweak Symmetry Breaking 

At some stage, we require that the electroweak symmetry is broken radiatively and use 
eq. m to determine the parameters /i^(MEwsB) and ^(Mewsb)- It is well known that 
the one-loop radiative corrections to the Higgs potential play a major role in determining 
the values of these two parameters, which at tree level are given in terms of the soft SUSY- 
breaking masses of the two Higgs doublet helds. We treat these corrections using the tadpole 
method. This means that we can still use eq. (CD) to determine /^^(Mewsb), one simply has 
to add one-loop tadpole corrections tu, td [22,25] 

m]j^ - tu/vu and ^ m]j^ - td/vd (39) 

We include the dominant third generation fermion/sfermion loops, as well as sub-dominant 
contributions from sfermions of the hrst two generations, gauge bosons, the MSSM Higgs 
bosons, charginos and neutralinos^, with the running parameters evaluated at Mewsb- Note 
that we also include the leading two-loop tadpole corrections, as calculated in Refs. [30-32], 
to be consistent in particular with the subsequent calculation of the physical neutral Higgs 
masses including the leading two-loop contributions. 

As far as the determination of /i^ and R/i is concerned, this is equivalent (at the considered 
loop level) to computing the full effective potential at scale Mewsb- Since |/r| and B affect the 
masses of some (s)particles appearing in these corrections, this gives a non-linear equation 
for |/i| (see Fig. 1), which is solved by a standard iteration algorithm until stability is reached 
and a consistent value of /i is obtained. From a practical point of view, this requires only 
three or four iterations for an accuracy of 0(10“^), if one starts from the values of |/r| and B 
as determined from minimization of the RG-improved tree-level potential at scale Mewsb 
and the procedure is extremely fast in CPU as compared to the (iterated) RGE calculation. 

At this stage, SuSpect includes a check on whether the complete scalar potential has 
charge and/or color breaking (CCB) minima which can be lower than the electroweak min¬ 
imum, or whether the tree-level scalar potential is unbounded from below (UFB). In the 
present version of the code, we consider only the following simple (tree-level) criteria [40] 

CCBl : Aj < 3 (ruj^-I--I--I-m^^). (40) 

UFBl ; > 2|R/i| at scale > M^-y^^gg (41) 

where / denotes any of the three fermion generations. Eq. iOD ensures that there is no deep 
CCB breaking minimum [due to very small Yukawa couplings] in some D-flat directions. 
One can either take this as a consistency necessary constraint on the MSSM parameters, 
or disregard it appealing to the fact that such minima are usually well separated from the 
electroweak minimum so that the latter can be reasonably stable at cosmological scales"^. 

^The contributions of the charginos and neutralinos can be rather sizable and are very important to 
minimize the scale dependence of the one-loop effective potential [24]. 

^But one would still lack for a compelling reason why the EW minimum is chosen in the first place. 
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For the third generation and in particular in the top sector, the CCB minimum is not much 
deeper than the electroweak minimum, since Yt is not very small, and not much separated 
from it. In this case one should apply eq. (P|) with some caution since tunneling effects can 
be important. On the other hand the “boundedness-from-below” condition of eq. (ED is 
actually an indication of possible dangerous non physical minima which could form when 
radiative corrections are included. At any rate, since both eqs. (HOD and ED are merely tree- 
level conditions, they should be checked at the highest energy scale. Note that all subsequent 
calculations (i.e. the sparticle and Higgs masses) are still performed, even if these conditions 
are not fulhlled, but a warning flag is given in the output hie. An upcoming version of 
SuSpect will have more sophisticated treatments, taking into account loop corrections [71] 
as well as the geometric conhgurations of the true minima [72,73] as will be discussed later. 

Finally, SuSpect provides appropriate (rejection) hags for any input choice in the pa¬ 
rameter space which lead to tachyonic pseudo-scalar Higgs boson or sfermion masses: 

NoTachyon: > 0 , > 0. (42) 

The electroweak symmetry breaking mechanism is assumed to be consistent when all these 
conditions are satished. Note however that for some input parameter choice such problems 
may occur temporarily i.e. before a sufficiently stable SUSY spectrum is obtained, and 
accordingly the SuSpect algorithm allows the calculation to be nevertheless performed until 
the very last iteration before eventually issuing a rejection hag. 

3.3 Calculation of the physical particle masses 

Once all the soft SUSY-breaking terms are obtained and eventually EWSB is radiatively 
realized [as should be the case in unihed models] one can then calculate all the physical 
particle masses. As already mentioned the whole procedure (namely, RGE -|- EWSB -|- 
spectrum calculation) is iterated a number of time until stability is reached (see the overall 
algorithm in Fig 1), in order to take into account realistic and stable particle masses in the 
“threshold” corrections using the expressions given in Ref. [25], to the gauge and Yukawa 
couplings at the scale Q = Mz- 

Our conventions for the mass matrices in the gaugino, sfermion and Higgs sectors will 
be specihed below. We basically follow the conventions of PBMZ with some important 
exceptions: {i) The fi parameter is dehned with the opposite sign (see below), {ii) The vevs 
are different by a factor a/2 and in our case v = 174.1 GeV. (Hi) The sfermion masses are 
dehned such that fi and /2 are, respectively, the lightest and the heaviest one. {iv) The 
matrices diagonalizing the chargino and neutralino mass matrices are taken to be real. 

For the calculation of the physical masses and the implementation of the radiative cor¬ 
rections, the various sectors of the MSSM are then treated as follows [with some details on 
the notation and conventions we use]: 
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3.3.1. The sfermion sector 

In the third generation sfermion sector [t,6, f], mixing between “left” and “right” current 
eigenstates is included [74], The radiatively corrected running DR fermion masses [essentially 
the Yukawa coupling times vevs] at scale Mewsb are employed in the sfermion mass matrices 
[this is important at large tan/?, where these corrections can be quite sizable]. As mentioned 
above, contrary to PBMZ, the masses are dehned such that and correspond to the 
mass of respectively, the lightest and the heaviest sfermion and therefore care should be made 
in interpreting the sfermion mixing angle 6^ as compared to PBMZ. [Note that a protection 
which prevents negative mass squared for third generation sfermions in the presence of too 
large mixing is provided.] The sfermion mass matrices are given by: 



+ (/| — e/s^)M| cos 2/5 + mf{Af — j^Vf) 

mf{Af — firf) cos2/3 + 


(43) 


where Af, fj, and rrif are respectively, the DR soft SUSY scalar masses, trilinear 

couplings, higgsino mass parameter and running fermion masses at the scale Mewsb and 
ri) = Tr = 1/rt = tan/5. These matrices are diagonalized by orthogonal matrices; the mixing 
angles 6f and the squark eigenstate masses are determined uniquely, at the relevant energy 
scales® by 


tan 26*/ = 


2 m/(A/ - /rr/) 
m| — m| + /rM| cos 2(3 

JL JR J 


(44) 


or 


tan 6*/ = 


m p — m‘3 + mA — m ^ — J?M| cos 2(3 

Jl _72_ JR _ JL J ^ __ 

2m,( A; - lit;) 


(45) 


and 
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™/ + - 


4 




y/(m3 — + ljM^cQs2jy-^ + — iir,)^ 


(46) 


Our convention for the mixing angles is 


A 

A 


cos 6*/ sin 6*/ 
— sin 6*/ cos 6*/ 



(47) 


where fi denote the mass eigenstates, and fL,R the chiral states. Equations (45, 46) ensure 
simultaneously that each mj, given in (44) is indeed the mass corresponding to the state 
fi and that /i is always the lightest state. In practice, the possible ambiguity, when Of is 

^Note that the sfermion masses and mixing angles are also calculated for consistency by SuSpect at the 
Mz scale as well since they enter the various radiative corrections to the gauge and Yukawa couplings. 
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determined from | arctan(tan26'/), is lifted in SuSpect by requiring that tan^j has always 
an opposite sign to A/ — /rr/, as can be seen from eq. (45)®. 

The radiative corrections to the sfermion masses are included according to Ref. [25], 
i.e. only the QCD corrections for the superpartners of light quarks [including the bottom 
squark] plus the leading electroweak corrections to the two top squarks; the small electroweak 
radiative corrections to the slepton masses [which according to PBMZ are at the level of one 
percent] have been neglected in the present version. 


3.3.2 The gaugino sector 


The 2x2 chargino and 4x4 neutralino mass matrices depend on the DR parameters 
Ml, M 2 , /i at the scale Mewsb and on tan/?. The chargino mass matrix given by: 


Me = 


M 2 \/ 2 Mw sin P 

\/ 2 Mw cos /5 n 


(48) 


is diagonalized by two real matrices U and V. The chargino masses are obtained analytically, 
with the convention that Xi is the lightest state. 

The neutralino mass matrix, in the M®) basis, has the form 


Ml 0 —MzSw cos (3 MzSw sin (3 

0 M 2 MzCw cos (3 —MzCw sin (3 

—MzSw cos P MzCw cos P 0 —/i 

MzSw sin P —MzCw sin P —jj, 0 


(49) 


It is diagonalized using analytical formulae [75] by a single matrix Z which is chosen to be 
real, leading to the fact that some (in general one) of the neutralino eigenvalues is negative. 
The physical masses are then the absolute values of these eigenvalues with some reordering 
such that the neutralinos X 123 A heavier with increasing subscript and Xi is the lightest 
neutralino. 

For the gluino, the running DR mass rrig at scale Mewsb is identihed with M 3 (Mewsb) 


mh®® — M3 (Mewsb) (50) 

The full one-loop QCD radiative corrections to the gluino mass are incorporated [20], while 
in the charginos/neutralinos case the radiative corrections to the masses are simply included 
in the gaugino and higgsino limits, which is a very good approximation [25]. 

®In other parts of the code where sfermion masses and mixing angles are used internally, for instance 
when entering in loop corrections, only the first equation in (44) is used to determine the mixing angles. In 
this case the consistency is guaranteed through the calculation of the eigenmasses directly from the mixing 
angles in the form (M| )ii cos^ 6f + (M|)22 sir? 6f ± (M?) 12 sin 20/, rather than from (45). Of course, the 

convention < m?~ is not necessarily valid here. 
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3.3.3. The Higgs sector 


In the latest version of SuSpect we have performed an up-to-date and more precise deter¬ 
mination of the physical Higgs masses including, in particular, the leading two-loop con¬ 
tributions [30-33] to the neutral Higgs masses and tadpoles. The running DR mass of the 
pseudoscalar Higgs boson at the scale Mewsb, Ma, is obtained from the soft SUSY-breaking 
Higgs mass terms evolved from RGE at the scale Mewsb and including the loop tadpole cor¬ 
rections [25] 


Ma{Mewsb) 


cos 2(3 


m 


H, 


U 

Vd 


m 


Hu 


-I- — I — m| -|- sin^ (3— + cos^ (3 — (51) 

Vd Vu 


This mass, together with the Z boson mass Mz at scale Mewsb, are then used as inputs in 
the CP-even Higgs boson 2 x 2 mass matrix Ais- Including the dominant contributions of 
the self-energies of the unrotated CP-even neutral Higgs helds and (as well as the 
tadpole contributions), this matrix reads at a given scale 


M| cos /32 +_Ml sin^ p - sn{q^) -\{Ml + M|) sin 2p - Si2{q^) 
-i(M| + Ml) sin 2p - si2(g2) M| sin^ p + Ml cos^ p - S 22 (g') 


One obtains the running CP-even Higgs boson masses in terms of the matrix elements Af k 

1 


Min = V 


2 L 




tfi 


Aff2)2 + 4(Af^ ^2 


'■ 12 ) 


(53) 


The mixing angle a which diagonalizes the matrix and rotates the helds into 

the physical CP-even Higgs boson helds h, H 


( H \ _ f cosa sino; \ \ 

\ h ) y — sin CK cos a y y j 


is given by 


sin 2 q; 


2 Aff 2 

Ml-Ml 


cos 2 q; 


Aif,-Aig 

Ml-Ml 


IT IT 
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(54) 


(55) 


The running charged Higgs boson mass at the EWSB scale is given by 


M^± = 


Ml + Mw 


• 2 ,3 

sm P — 
Vd 


2 n 

COS p — 

V„ 


(56) 


The pole masses of all the Higgs bosons are then obtained by including the self-energy 
corrections evaluated at the masses of the Higgs bosons themselves. 

In the evaluation of the radiative corrections in the MSSM Higgs sector which are known 
to be very important [27], we have made several options available: 

(f) Approximate one-loop and two-loop contributions to the self-energies (and tadpole) 
corrections Sij in the mass matrix Ai^■ These expressions, given in Ref. [29], provide a rather 
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good approximation (at the few percent level) for the masses and Mh and the angle a 
at least for a reasonable range and not too extreme values of input parameters. Since it 
makes the program running faster, this approximation may be convenient in cases where a 
very precise determination of the Higgs masses may not be mandatory. 

{a) A full one-loop calculation of the Higgs masses, following essentially the expressions 
of Ref. [25]. In this case, the tadpoles are consistently evaluated at the one-loop level only. 

iiii) A full one-loop plus leading two-loop calculation controlled by the third generation 
Yukawa couplings and the strong gauge coupling, derived in Refs. [30-33]. Note that the 
calculation is entirely performed in the DR scheme, therefore different from the on-shell 
scheme calculation performed in FeynHiggsFast of Heinemeyer, Hollik and Weiglein [28,29]. 
The comparison of results from {ii) and {Hi) can be useful in order to see the effects from 
two-loop contributions. For a recent detailed analysis of the Higgs masses in various models 
and a comparison between SuSpect and other codes, we refer to Ref. [34]. 

3.4 Theoretical and experimental constraints on the spectra 

Once the SUSY and Higgs spectrum is calculated, one can check that some theoretical and 
experimental requirements are fulhlled. Examples of theory requirements are for instance, 
the absence of charge and color breaking (CCB) minima and that the potential is not un¬ 
bounded from below (UFB), the absence of too much hne-tuning (FT) in the determination 
of the masses of the Z boson from EWSB as well as in the determination of the top quark 
mass. For experimental requirements on the spectrum, one can demand for instance that it 
does not lead to large radiative corrections to the precisely measured electroweak parameters 
or too large values for the anomalous magnetic moment of the muon and the branching ratio 
of the radiative decay of the fe-quark. SuSpect provides such tests in the form of warning 
or error flags in the output hie. 

3.4.1 CCB and UFB 

As explained previously, in Suspect, the EWSB conditions are consistently implemented 
by iteration on the parameters p and B and the occurrence of a local minimum is checked 
numerically. In the same time one needs to check for the non existence of deep CCB minima 
or UFB directions. Avoiding such cases may put strong constraints on the model and we 
mentioned in section 3.2.3 that we have already implemented two simple CCB and UFB 
conditions [40] and the program gives a warning hag when they are not satished. 

In a next version of the code, to be released in a near future, we will address the question 
of CCB minima and UFB directions in the most complete possible way, given the present 
state of the art. Three complementary features should be considered in relation to the CCB 
minima: {i) the directions in the space of scalar helds along which such minima can de¬ 
velop {ii) whether they are lower than the EWSB minimum {Hi) whether the EWSB (then 
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false) vacuum can still be sufficiently stable. In Ref. [71] a systematic study of point (i) 
has been carried out considering subspaces involving the helds Hu, Qu, ur {Hd and possibly 
L). However, the identified D-flat directions contain the true minima only in the case of 
universal scalar soft masses at the low energy relevant scales, otherwise they constitute only 
sufEcient conditions for the occurrence of CCB minima. While such directions provide very 
good approximations for the hrst two generations, special attention should be paid to the 
third generation sector as was stressed in Ref. [72]. This is relevant in particular to codes 
like SuSpect where various SUSY model assumptions can be considered, including non- 
universality. Furthermore, the check of point {ii) as done in Ref. [71] involves a numerical 
scan over held values. Actually there are cases where held-independent conditions can be 
obtained even in the case of 5-held directions Hu, Q,ur, Hd, L, leading to faster algorithms; 
see Ref. [72] and unpublished study. We will thus optimize in SuSpect the various avail¬ 
able complementary approaches. Point {Hi) has also its importance as it can increase the 
phenomenologically allowed regions of the MSSM parameter space. Some simple criteria 
will be encoded, following for instance Ref. [73]. Finally, the UFB directions as identihed 
in Ref. [71], in particular UFB-3, lead to very strong constraints. Nonetheless, there is still 
room for some improvements by optimizing the criterion of “deepest direction”, leading in 
some cases to even stronger constraints which will be also implemented in SuSpect. 


3.4.2 Fine-tuning 


One of the main motivations for low energy SUSY is that it solves technically the hierarchy 
and naturalness problems. However, since the Z boson mass is determined by the soft SUSY- 
breaking masses Mjj^, and the parameter fj? as can be seen from eq. (11), naturalness 
requires that there are no large cancellations when these parameters are expressed in terms of 
the fundamental parameters of the model [for instance mo, mi/ 2 , 1^, B in mSUGRA], otherwise 
hne-tuning is re-introduced [19,77,78]. A similar problem occurs in the case of the top quark 
mass, since it is related to the top Yukawa coupling and tan /?. Various criteria for quantifying 
the degree of fine-tuning in the determination of Mz and m^ have been proposed and some 
subjectivity is involved in the statement of how much hne tuning can be allowed. Therefore, 
in our case, we simply evaluate the sensitivity coefficients for M| and m^ with respect to a 
given parameter a [78] 


5Ml 

Smt 

mt 


C{Ml,a)- 

a 

C{mt,a) — 
a 


(57) 


and leave to the user the decision of whether the amount of hne-tuning [large values of the C 
coefficients] is bearable or not. We evaluate only the hne-tuning with respect to variations 
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of the parameters jj? and Bjj,, for which the coefficients take the simple form: 


FTIMZ : 

. ^^,^2 ,2^ 2p2 ( 4 p{ml-ml)' 


FTIMZ : 



FTIMT : 



FTIMT : 


(58) 


with tfs = (tan^/3 + l)/(tan^/3 — 1). Further hne-tuning tests can be made, in particular 
with respect to the f, b Yukawa couplings and are planned to be included in future versions. 

3.4.3 Electroweak precision measurements 

Loops involving Higgs and SUSY particles can contribute to electroweak observables which 
have been precisely measured at LEP, SLC and the Tevatron. In particular, the radiative 
corrections to the self-energies of the W and Z bosons, Ilww and Hzz, might be sizable if 
there is a large mass splitting between some particles belonging to the same SU(2) doublet; 
this can generate a contribution which grows as the mass squared of the heaviest particle. 
The dominant contributions to the electroweak observables, in particular the W boson mass 
and the effective mixing angle enter via a deviation from unity of the p parameter [79], 
which measures the relative strength of the neutral to charged current processes at zero 
momentum transfer, i.e. the breaking of the global custodial SU(2) symmetry: 

p = (1 — Ap)“^ ; Ap = n^2(0)/M| — nww(0)/M-^ (59) 

Most of the MSSM contributions to the p parameter are small, Ap < 10“^ [80]. In the case 
of the Higgs bosons, the contributions are logarithmic, ~ Q;log(M/j/M^), and are similar 
to those of the SM Higgs boson [and identical in the decoupling limit]. The chargino and 
neutralino contributions are small because the only terms in the mass matrices which could 
break the custodial SU(2) symmetry are proportional to Mw- Since in general, first/second 
generation sfermions are almost degenerate in mass, they also give very small contributions to 
Ap. Therefore, only the third generation sfermion sector can generate sizable corrections to 
the p parameter, because of left-right mixing and [in case of the stop] the SUSY contribution 
oc m'j leads to a potentially large splitting between the sfermion masses. 

We have thus calculated Ap in the MSSM, taking into account only the contributions of 
the third generation sfermions. We include full mixing and in the case of the stop/sbottom 
doublet, also the two-loop QCD corrections due to gluon exchange and the correction due 
to gluino exchange in the heavy gluino limit, which can increase the contribution by 30% 
or so [81]. One may then require the SUSY contribution not to exceed e.g. two standard 



deviations from the SM expectation [82]: Ap(SUSY) <2-10 The precise value of Ap as 
calculated by SuSpect for any choice of parameters is written in the output hies. 

3.4.4 The muon (g-2) 

The muon {g — 2) anomalous magnetic moment has been very precisely measured and thus 
generally can set strong constraints on the additional contribution from SUSY particles. The 
latest experimental value is [83]: 

{g^ - 2) = a;fP = (11 659 208 ± 6) 10“^°, (60) 

The theoretical value predicted in the SM, including the QED, QCD and electroweak correc¬ 
tions, has some uncertainties from the determination of the hadronic vacuum polarization 
contributions. In particular, at the time of this writing, the SM theoretical prediction de¬ 
pends quite crucially on whether one takes into account the combined measurement from 
e^e“ and r decays data, or from e^e~ data only, for the calculation of the hadronic vacuum 
polarization via a dispersion relation [84]. 

The contribution of SUSY particles to (p^ —2) [85,86] is mainly due to neutralino-smuon 
and chargino-sneutrino loops [if no havor violation is present as is the case here]. In typical 
models (such as mSUGRA), the contribution of chargino-sneutrino loops usually dominates. 
This is also true for conhgurations where all superpartners are almost degenerate in mass. 
In the latter case the contribution of loops can be approximated by [m is the mass 

of the heaviest particle per GeV]: ~ 10“^ x (tan/5/m^), to be compared with the 

contribution of loops, Aa^°^ ~ —10“® x [iaxi (3/fh?). These contributions are large for 
large values of tan (3 and small values of the scalar and gaugino masses. Moreover, in some 
other parts of the parameter space and for relatively heavy SUSY particles, the dominant 
SUSY effect follows the sign of the p parameter. 

We have included a routine which calculates the full one-loop contributions of chargino- 
sneutrino and neutralino-smuon loops in the MSSM, using the analytical expressions given 
in Ref. [86]. In this case, the full mixing in the smuon sector is of course included [this is the 
only place where the parameter plays a role in the code]. We also took into account in 
this evaluation the leading two-loop QED correction [87], which essentially reduce by about 
~ 7% the one-loop contribution. 

In view of the above mentioned present uncertainties of the SM predictions, we refrain 
from providing a particular interpretation of the comparison between the experimental above 
value and one of the SM theoretical predictions, in terms of e.g. standard deviations from 
SM. Rather, the supersymmetric contributions are simply written in the SuSpect output 
hies for any choice of input parameters and models, and can be directly compared with the 
(eventually updated) preferred range of the quantity 
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3.4.5 The radiative decay b ^ s'y 

Another observable where SUSY particle contributions might be large is the radiative flavor 
changing decay b ^ sj [88,89]. In the SM, this decay is mediated by loops containing charge 
2/3 quarks and hh-bosons but in SUSY theories, additional contributions come from loops 
involving charginos and stops or top quarks and charged Higgs bosons [contributions from 
loops involving gluinos or neutralinos are very small [88] in the models considered here]. 
Since SM and SUSY contributions appear at the same order of perturbation theory, the 
measurement of the inclusive B Xg'y decay branching ratio [41] is a very powerful tool 
for constraining the SUSY parameter space. 

In Refs. [89,90], the next-to-leading order QCD corrections to the decay rate in the 
MSSM have been calculated and a Fortran code has been provided. It gives the most up- 
to-date determination of BR(6 —sq) where all known perturbative and non-perturbative 
effects are implemented, including all the possibly large contributions which can occur at 
NLO, such as terms oc tan/? and/or terms containing logarithms of Mewsb/^w- We have 
interfaced this routine with our code^. Besides the fermion and gauge boson masses and 
the gauge couplings that we have as inputs, we will use the values of the other SM input 
parameters required for the calculation of the rate given in Ref. [91], except for the cut-off 
on the photon energy, > {1 — S)mi,/2 in the bremsstrahlung process b s'yg, which we 
fix to (5 = 0.9 as in Ref. [90]. Then, again, the theoretical value of BR(6 —sq) as calculated 
by the above mentioned routine is written in the output files, and can be compared by the 
user with the latest available data [41]. 


4. Running SuSpect 

4.1 Basic facts about SuSpect 

The program Suspect is composed of several files and routines: 

i) The input files suspect2. in (standard format) or alternatively suspect2_lha. in (the 
SLHA format): here one can select the model to be investigated, the accuracy of the algo¬ 
rithm, the input data (SM fermion masses and gauge couplings). Some reasonable default 
values are set in the example which is provided. One would then simply select the SUSY 
model (pMSSM, mSUGRA, GMSB and AMSB), choose the corresponding input parameters 
and possibly make a few choices concerning the physical calculation (such as enforcing or 
not unification of the gauge couplings, changing the scale at which EWSB occurs, including 
or not RG to the masses, and choosing the routine calculating the Higgs masses). The list 
of possible choices is given in the next subsection. 

ii) The program suspect2_call. f: this is an example of a routine which calls the main 
subroutine suspect2.f. This program is necessary to run SuSpect since it defines the 

^We thank Paolo Gambino for providing us with his code and for his help in interfacing it with ours. 
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primary algorithm control input parameters needed by the latter. In particular, there is an 
important parameter (INPUT) which controls the form of input/output: depending on its 
value it allows to select the input and output either in the standard (suspect2. in) or in the 
SLHA (suspect2_lha. in) format, or alternatively to bypass the reading of any input hies, 
in which case all the parameters and choices are to be dehned by the user in this calling 
routine. This last choice is particularly useful for interfacing SuSpect with other routines 
and/or for scans of the parameter space. 

Hi) The main routine suspect2.f: here all the calculation of the spectrum is performed, 
once the input is supplied by suspect2.in, suspect2_lha. in or suspect2_call.f. This 
routine is self-contained, except for the determination of the Higgs boson masses where 
it needs to call one additional routine for the two-loop contributions that we also supply: 
twoloophiggs. f (originally provided by Pietro Slavich). Note that we provide also as a 
separate hie, needed in the compilation: bsg.f, performing the calculation of the b ^ s'j 
branching ratio (originally provided by Paolo Gambino). 

iv) output hies suspect2.out and suspect2_lha.out: these hies are generated by de¬ 
fault [they can be switched oh by an appropriate value of the control parameter INPUT 
in suspect2_call. f] at each run of the program and gives the results for the output soft 
SUSY-breaking parameters [when they are calculated] and the masses and mixing angles 
of the Higgs and SUSY particles. Some warnings and comments are also given when the 
obtained spectrum is problematic as will be discussed in the next subsection. 

The routine suspect2.f consists of about 10.000 lines of code and takes about 350 Ko 
of memory, while the input and the calling routines have only a few hundred lines (most of 
them being comments). The accompanying separate routine for the calculation of the two- 
loop radiative correction to the Higgs masses has about 7000 lines of code. The complete 
executable code takes about 0.9 Mo space. 

The Fortran hies have to be compiled altogether and, running for instance on a PC 
using GNU Fortran, the (minimal options) compilation and link commands are: 
g77 -c suspect2_call.f suspect2.f twoloophiggs.f bsg.f 
g77 -o suspect suspect2_call.o suspect2.o twoloophiggs.o bsg.o 
suspect 

no other compilation option is in principle needed, though some users might hnd it useful 
to try other standard compilation options. The running time for a typical model point, for 
instance an mSUGRA point, is about 1 second on a PC with a 1 GHz processor. 

4.2 The main routine and the control parameters 

In this subsection, while we will refrain from exhibiting the input and output hies (some 
examples can be found on the web page of the program), we mention a few important 
features about the main routine and the control parameters that it is uses. 
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SUBROUTINE SUSPECT2(iknowl,input,ichoice,errmess) 


is the main routine of the program, to be used as it is or to be called by any other routine 
(such as suspect2_call. f, as will be discussed below). It has the following four basic input 
control parameters: 

IKNOWL; which sets the degree of control on the various parts of the algorithm. It has 
two possible values: 

- IKN0WL=0; blind use of the program, i.e. no control on any “algorithmic” parameter. 
Reasonable default values are set for the control parameters and the program gives 
just the results from the physical input. 

- IKN0WL=1; in addition warning/error messages are collected in the suspect2.out hie 
(this is the recommended choice in general). 

INPUT; is for the physical input setting and can take the following values: 


- INPUT=0; the model and option parameters ichoice(l)-( 10) as well as the val¬ 
ues of the input parameters are read off from the hie suspect2.in in the original 
SuSpect format. The output is written in both suspect2.out and SLHA format 
suspect2_Iha.out hies. 

- INPUT=1: the user dehnes all the relevant input choices and parameters within his 
calling program, i.e. any reading of input hle(s) is bypassed. The required list of 
parameters to be dehned (with consistent names etc), can be found in the commons 
given below. In addition, examples of input parameter setting for this option, with 
all the necessary input parameters (appearing in common blocks), are given for vari¬ 
ous models in the calling program suspect2_call. f. The output is written in both 
suspect2.out and SLHA format suspect2_Iha.out hies. 

- INPUT=2: the model and option parameters ichoice(l)-( 10) as well as the val¬ 
ues of the physical input parameters are read oh from the hie suspect2_lha. in in 
the SLHA format. The output is written in both suspect2.out and SLHA format 
suspect2_lha.out hies. 

~ INPUT=11; same as INPUT=1, but with no output hle(s) written (note this option may 
be more convenient e.g. for scans of the MSSM parameter space). 


ICHOICE; initializes the various options for the models to be considered, the degree of 
accuracy to be required, the features to be included, etc. There are 11 possible choices at 
present and the options are described in detail in the input hies suspect2. in or suspect2_lha. in; 


- ICHOICE(l) 
~ ICH0ICE(2) 

- ICH0ICE(3) 
~ ICH0ICE(4) 
~ ICH0ICE(5) 

- ICH0ICE(6) 


Choice of the model to be considered. 

For the perturbative order (1 or 2 loop) of the RGBs. 
To impose or not the GUT scale. 

For the accuracy of the RGBs. 

To impose or not the radiative BWSB. 

To chose diherent input in general MSSM. 
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~ ICH0ICE(7): For the radiative corrections to the (s)particles masses. 

“ ICH0ICE(8); To set the value of the EWSB scale. 

- ICHOICEO); For the hnal accuracy in the spectrum calculation (two possibilities). 

- ICHOICE(IO); For the approximation in calculating the Higgs boson masses (approxi¬ 
mate, one-loop, or two-loop). 


ERRMESS; which 
cally written in the 

- ERRMESS(i)= 

- ERRMESS(1)=- 

- ERRMESS(2)=- 
~ ERRMESS(3)=- 
~ ERRMESS(4)=- 
~ ERRMESS(5)=- 

- ERRMESS(6)=- 

- ERRMESS(7)=- 
~ ERRMESS(8)=- 
~ ERRMESS(9)=- 

- ERRMESS(10)= 


provides a useful set of warning/error message flags, that are automati- 
output hies suspect2.out or suspect2_lha.out; 

0; Everything is hne, 

T; tachyonic 3rd generation sfermion from RGE, 

T: tachyonic lst/2d generation sfermion from RGE, 

T: tachyonic A boson (maybe temporary: see hnal mass), 

T; tachyonic 3rd generation sfermion from mixing, 

T: /i(MGUT) guess inconsistent, 

T; non-convergent fi from EWSB, 

T: EWSB maybe inconsistent (but RG-improved only check), 

T: ^iggs uiaybe UFB or GGB (but RG-improved only check), 

T: Higgs boson masses are NaN, 

=-l; RGE problems (non-perturbative and/or Landau poles). 


4.3. Calculations with SuSpect 

4.3.1 Comparison with other codes 

Our results for some representative points of the MSSM parameter space have been carefully 
cross-checked against other existing codes. Most of the earlier comparisons [94] have been 
performed in the context of mSUGRA type models but, more recently, comparisons in other 
models have been performed. We obtain, in general, a very good agreement, at the percent 
level and often better [95], with the codes SoftSusy and SPheno®. The most sophisticated 
parameter to obtain in this context is the lightest Higgs boson mass, since it incorporates all 
possible ingredients: the RGE’s for the evaluation of Mh„ and Mh^, the effective potential 
and the EWSB for the determination of Ma and the tadpoles, the radiative corrections to the 
Higgs sector which involve also the two-loop corrections, etc. Recently, a detailed comparison 
of SuSpect with SoftSusy and SPheno (and also FeynHiggs) for the Higgs boson masses in 
many scenarios (mSUGRA, GMSB, AMSB, general MSSM) was performed in Ref. [34]. The 
value that we obtain for for instance, is in a very good agreement, the difference being 
less than half a GeV, while the agreement for Ma and fi is at the level of 1%; see Ref. [34] 
for details. The agreement with the program ISASUGRA is not as good, in particular for what 
concerns the Higgs sector, since different approximations have been used in the two cases. 

®We thank Ben Allanach, Sabine Kraml, Werner Porod and Pietro Slavich for their gracious help in 
performing these detailed comparisons of the programs. 
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4.3.2 Interface with other programs 

In the way it is written, SuSpect can be easily interfaced with other programs or Monte- 
Carlo event generators®. Private or official versions of some programs exist in which an 
interface with SuSpect has been made; we give a short list of them^®: 

• micrOMEGAs [96]: for the automatic (analytical and then numerical) calculation of the 
cosmological relic density of the lightest neutralinos, including all possible channels. 

• DARKSUSY [97]: also for the calculation of the relic density of the lightest neutralinos 
[including co-annihilation effects] and their direct and indirect detection rates as well 
as other astrophysical features. [The program has its own calculation of the SUSY 
spectrum, but it is rather approximate]^^ . 

• HDECAY [76]: for the calculation of the decay branching ratios and total decay widths of 
the SM and MSSM Higgs bosons [in fact some routines, in particular those for the QCD 
running and for the interface with the routines calculating the Higgs boson masses, are 
borrowed from there]. 

• SDECAY [100]: for the calculation of the decay widths and branching ratios of SUSY 
particles including higher order [three-body decays for gauginos and stops, four-body 
decays for the lightest stop and QCD corrections to the two-body decays of squarks 
and gluinos], which will appear soon. 

• SUSYGEN [101]: a Monte-Carlo event generator for Higgs and SUSY particle production 
in the MSSM [mainly in e’''e“ collisions but some processes in ep and pp collisions are 
implemented]. The program is also interfaced with HDECAY. 

• SFITTER [102] which determines the weak scale MSSM parameters from measurements 
performed at hadron and colliders from a combined £t/grid taking into account 
the correlations. 

As discussed already, we also interfaced SuSpect with the code bsg.f which calculate 
the rate of the b ^ s'j decay at next-to-leading order [90]. An interface with the Monte- 
Carlo event generators such as PYTHIA [103] and HERWIG [104] can be made throgh the SLHA 
format. 

®To make this interfacing easier, in addition to including the input and output files in the now standard 
SLHA format, we also have provided a set of obvious commons for the input and output parameters needed 
or calculated by SuSpect and named all commons, subroutines and functions used with a prefix SU_, in order 
not to be in conflict with those used by other programs. 

^°Note that we have also interfaced the program with a private code written by Manuel Drees calculating 
the cosmological relic density of the lightest neutralinos for the complete analysis of the mSUGRA parameter 
space performed in Ref. [92]. 

^^Note that a private Suspect/DARKSUSY interface has been already used for prediction studies of indirect 
LSP detection [98,99]. 
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4.3.3 Web information and maintenance 


A web page devoted to the SuSpect program can be found at the http address: 

http: //www. Ipta. univ-montp2. fr/~ kneur/Suspect 

It contains all the information that one needs on the program: 

- Short explanations of the code and how to run it. 

- The complete “users manual” can be obtained in post-script or PDF form. 

- A list of important changes/corrected bugs in the code. 

- A mailing list to which one can subscribe to be automatically advised about future 
SuSpect updates or eventual corrections. 

One can also download directly the various hies of the program: 

- suspect2.in, suspect2_lha. in: the input hies in standard or SLHA format. 

- suspect2_call. f: the calling program sample. 

- suspect2.f: the main routine of the program. 

- twoloophiggs . f: the additional Higgs routine needed. 

- SuSpect2.tar.gz: all needed routines in standard compressed format for the latest 
version. 

- SuSpect2_01d.uu: the routines for the previous versions of the program. 

Note that the possibility to use directly the program [as well as the other RGE codes] 
interactively on the web has been set by Sabine Kraml [95]. 

5. Conclusion 

We have presented the version 2.3 of the Fortran code SuSpect which calculates the 
Supersymmetric and Higgs particle spectrum in the MSSM. The calculation can be performed 
in constrained models with universal boundary conditions at high scales such as the gravity 
(mSUGRA), anomaly (AMSB) or gauge (GMSB) mediated breaking models, but also in 
the non-universal or unconstrained MSSM case, with up to 22 free input parameters which 
can be set either at the electroweak symmetry breaking scale or obtained from boundary 
conditions on some common parameters at a high-energy scale. 

A particular care has been taken to treat all the mandatory features which are needed to 
describe accurately these various scenarii: the renormalization group evolution of parameters 
between low and high energy scales, the consistent implementation of radiative electroweak 
symmetry breaking and the calculation of the physical masses of the Higgs bosons and su¬ 
persymmetric particles taking into account all dominant radiative corrections. The program 
provides several options [for accuracy, scale choice, etc...] to deal with these aspects. 

The program can check the fulhllment of theoretical constraints, such as the absence of 
tachyonic particles and improper lightest SUSY particle, the absence of non desired charge 
and color breaking as well as unbounded from below minima and a large hne-tuning in the 
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electroweak symmetry breaking conditions. A verification that the obtained spectrum is in 
agreement with high precision measurements such as the p parameter, the muon g — 2 and 
the radiative b ^ s'j decay, can also be performed. 

The program has a high degree of flexibility in the choice of the model and/or the input 
parameters and an adequate level of approximation at different stages. It is rather precise and 
quite reliable [since it has been compared with several other similar existing codes], relatively 
fast to allow for rapid comprehensive scans of the parameter space and simple enough to 
be linked with other programs dealing with MSSM particle properties or with Monte-Carlo 
event generators. We have also provided a simpler way to run the code interactively on the 
web in the constrained models. 

Several upgrades, which include the possibility to analyze additional theoretical models 
and to further extend the already existing interfaces with programs for (s)particle decay 
branching ratios and Dark Matter calculations and with some Monte-Carlo event generators 
to simulate the production properties, are planned. 
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